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ABSTRACT 

We present period change rates (dP/dt) for 42 RR Lyrae variables in the globu- 
lar cluster IC4499. Despite clear evidence of these period increases or decreases, the 
observed period change rates are an order of magnitude larger than predicted from 
theoretical models of this cluster. We find there is a preference for increasing periods, 
a phenomenon observed in most RR Lyrae stars in Milky Way globular clusters. The 
period-change rates as a function of position in the period-amplitude plane are used to 
examine possible evolutionary effects in Ool clusters, OoII clusters, field RR Lyrae stars 
and the mixed-population cluster lo Centauri. It is found that there is no correlation 
between the period change rate and the typical definition of Oosterhoff groups. If the 
RR Lyrae period changes correspond with evolutionary effects, this would be in contrast 
to the hypothesis that RR Lyrae variables in OoII systems are evolved HB stars that 
spent their ZAHB phase on the blue side of the instability strip. This may suggest that 
age may not be the primary explanation for the Oosterhoff types. 

Subject headings: surveys — stars: abundances, distances. Population II — Galaxy: 
center 



1. Introduction 



Within our Galaxy, most globular clusters that contain RR Lyrae variable stars can be divided 
into two distinct groups: in Oosterhoff I clusters the fundamental-mode RR Lyrae variables have an 
average period near 0.55 days, while in Oosterhoff II clusters the average fundamental-mode period 
is close to 0.65 days. The two types of cluster also differ in the relative numbers of first overtone 
pulsators (RRl-type variables) versus fundamental mode pulsators (RRO-type variables), a nd in 
their period-amplitude (PA) relations (jLee Carnevlll999bl : IPritzl et al.ll20od : ICatelanlboog ). 



The Oosterhoff dichotomy has been a source of great interest for decades in modern stellar 
astrophysics, especially after it was shown that our Galaxy is the only known example to exhibit 
such a dichotomy. The RR Lyraes in neighboring dwarf galaxies have different values of < Prro > 



Based in part on observations made with the European Southern Observatory telescopes obtained from the 
ESO/ST-ECF Science Archive Facility. 
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than either of the Ool and OoII groups {e.g., IPairOra et alj|2003l : ICatelaiJl2004l . [20091; ICreco et al 
2003 ) . An understanding of the Oosterhoff effect is essential to any understanding of t he formation 



and evolution of the globu lar cluster system and the formation history of the Galaxy (jEggen et al, 
19621 : ISearle k ZinnlflQTsI ^ . 



The leading theoretical explanation for this phenomenon involves a dichotomy in the "transi- 
tion period" between th e RRO and RRl varia,bles w hich reflects a difference in effective temperature 
at the transition point ([van Albada &: Bakeiill97lh . This "hysteresis effect" would cause a delay in 
the mode switching (from e.g., RRO to RRl) and would occur at different temperatures, depending 
on the direction of evolution. According to this scenario, for the Ool, intermediate metallicity 
clusters, most of the horizontal-branch stars begin their lives in or near the instability strip. The 
RR Lyraes become hotter and their periods become shorter as they evolve and eventually the RRO 
stars switch their mode to become RRl-type variables. For the OoII, more metal-poor clusters, the 
HB stars begin their lives on the hotter (blue) side of the instability strip and subsequently evolve 
to lower temperatures and longer periods as they traverse the instability strip at luminosities ap- 
preciably higher than the zero-age horizontal branch. Thus, a luminosity (mean density) difference 
contributes to t he effect in addition to the direction of evolution as per van Albada and Baker. 
As an example, iBono et al.l (j 199 71 ) show evolutionary tracks of HB stars with the predicted insta- 
bility strip and point to where the transition is expected to take place according to the hysteresis 
hypothesis. 

This explanation also has the advantage that it would help explain the larger fraction of RRl 
variables in OoII clusters. The metallicity of the cluster is the "first" parameter that establishes 
where on the horizontal branch the average star settles after core helium burning begins. But the 
exact ZAHB location of each individual star is governed by the efficiency of the mass-loss along 
the RGB. The more metal-poor clusters have bluer horizontal branches and would become OoII 
clusters and vice versa for the intermediate metallicity Ool clusters. However, this explanation 
does not explain why no dichotomy is seen in the globular clusters belonging to dwarf spheroidal 
galaxies unless the range in metallicity of the old population in each galaxy is very small. 

Although it is common to categorize the GCs as Ool and OoII clusters, the Oosterhoff effect 
can manifest itself on a star-by-star basis, rather than only for the mean values of a given cluster 
(Sandage, Katem Sz Sandage 1981). The period-amplitude relations for RR Lyrae variables in Ool 
and OoII globular clusters differ substantially, and the position of an RR Lyrae in the PA plane is 
often used as a diagnostic to determine whether an RR Lyrae is Ool- or OoII-type. For example, 
it is not uncommon to find a few RR Lyrae variables in Ool globular clusters that occupy the OoII 
PA position and these are believed to be more evolved and intrinsically brighter than the majority 
of the Ool stars {e.g., Cacciari, Corwin & Carney 2005, Sandage et al. 1981). If evolution is the 
explanation for the Oosterhoff groups, then this shift of the stars' positions in the PA plane by 
Ool and OoII RR Lyrae stars would be explained not by a simple difference in abundances, (see 
Clement & Shelton 1999), but by evolutionary effects. However, it was shown that both the "young" 
and "old" halo GCs (using the Mackey &: van den Bergh 2005 classification scheme) present the 
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Oosterhoff dichotomy (|Catelanll2009l ). in contrast to the notion that OoII clusters as a group are 
more evolved and hence older than Ool clusters. When using the (smaller-sample) Marfn-Franch 
et al. 2009 relative age classification scheme, this is not as clear. 

Another explanation for the Oosterhoff groups is that helium abundance is affec ting the RR 



Lyrae vari ables in the Ool and Poll globular clusters. It has recently been shown by iBusso et al 



(|2007l ) and lCaloi &: D'Antonal (|2007l ) that a helium enriched population could explain the HBs in 
the peculiar, metal-rich globular clusters NGC 6388 and NGC 6441. These two metal-rich GCs 
([Fe/H]~— 0.6 dex) have complicated the Oosterhoff groups since they have been found to contain 
RR Lyrae variables with unusually long periods. In terms of the evolutionary explanation for the 
Oosterhoff groups, th ese would be clearly classified as extreme Oosterhoff type II clusters. However, 
it has been shown by lPritzl et al.l (120001 ) that neither their HB morphology in the color-magnitude 
diagrams nor theoretical models by Sweigart & Catelan (1998) indicate that the RR Lyraes are 
in a more advanced evolutionary stage. Instead, the presence of multiple stellar populations with 
different initial He contents could account for the peculiar morphology of the HB in both clusters 
(jBusso et al.l 120071 ). Further, it is seen that clusters with different metallicities exhibit similar HB 
morphologies, while clusters of the same metal abundance can have rather different HB morpholo- 
gies. This is the "second parameter" puzzle, which also may be explained (at least in part) by 
differences in helium content (see e.g., Gratton et al. 2010). 

As a by product of understanding the Oosterhoff effect, vital information on the second pa- 
rameter phenomenon may be obtained. We will here investigate the Oosterhoff dichotomy by 
comparing the evolutionary status of Ool and OoII RR Lyrae stars, where the evolutionary status 
is determined from the RR Lyrae period change rates. In addition we report on a pilot study to 
formulate Period-Color (PC) and Amplitude-Color (AC) relations as a function of pulsation phase 
and Oosterhoff type. These resulting relations could permit a more detailed comparison of RR 
Lyrae properties with Oosterhoff type. 

New period change rates are obtained for the RR Lyrae variables in IC 4499. IC 4499 is consid- 
ered an Ool cluster, with a mean RR Lyrae period of < Prro >=0.579 d and < Prri >=0.347d, 
and a ratio of the number of RRl stars to the total number of RR Lyrae of Nrri/Nrr=0.20. The 
IC4499 period change rates are compared to the period change rates in the Ool clusters M5, M3 
and NGC 7006, the OoII clusters M15, M2, NGC 5466, M22 and NGC 5053, the field RRO Lyrae 
variables and the mixed population cluster cj Cen. Period change rates give insight on the evolu- 
tion of RR Lyrae variables. The periods of the RR Lyrae variables should be increasing if the star 
evolves from b lue to red in the HR diagram and decreasing if the star evolves from red to blue 
(jSandagd 119571 ). Hence, if evolution is responsible for the difference between Ool and OoII- type 
RR Lyrae variables, those stars with Ool properties would evolve from red to blue and should have 
decreasing periods. Similarly, those RR Lyrae stars with OoII properties should ha ve periods tha t 
increase. However, RR Lyrae period changes are not fully explained by evolution (jRosindll973l ). 
and it is frequently seen that observed peri od changes are an order o f magnitude larger than would 
be expected from HB evolutionary models (jSweigart &: Renzinilll979l ). 
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By comparing where the increasing and decreasing RR Lyrae period change rates fall on the 
PA diagram, we can explore whether a relation exists between Oosterhoff type and period change 
rate (evolutionary state?). Although mean RR Lyrae period change rates have been calculated for 
different globular clusters of different Oosterhoff groups, the individual RR Lyrae period change 
rates (and hence direction of evolution) as a function of position on the PA plane has not been 
explored. 

This paper first addresses the observational data for the RR Lyrae variables in IC4499 in 
§2. The determination of their period change rates follows in §3 and the presentation of these 
period change rates is given in §4. Next, in §5, the dereddened ( V-R) colors and minimum light 
temperatures of the IC4499 (Ool) RRO variables are compared to the RRO in M15 (OoII). In §6, 
implications regarding the period change rates of RRO Lyrae variables in the different Oosterhoff 
GCs are discussed. A summary is presented in §7. 



2. Data and Photometry 

The present photometric data reduction for IC4499 comes from analysis of 1,338 individual 
CCD images in the direction of IC 4499 from 23 observing runs on five telescopes (Magellan/Baade 
6.5m, CTIO Blanco 4m, ESO NTT 3.6m, CTIO 1.5m, and CTIO 0.9m). These data have been 
placed by one of the authors (PBS) on the photometric system of Landolt (1992; see Stetson 2000, 
2005), and span a range of slightly more than 22 years (1987 March 2 to 2009 April 3). The various 
cameras used had different physical sizes and different image scales, therefore slightly different 
areas around IC 4499 were acquired and any given star may not appear in all images. The average 
number of independent observations for a star was 165 in B, 244 in V, 36 in R and 90 in /. More 
details on the data sets can be found in our paper discussing the IC 4499 color magnitude diagram 
(Walker et al, in preparation). 

Our analysis of the RR Lyrae light curves also makes use of the original photographic B 
data from Clement, Dicken s & Bingham(1979). The V and R photometry for Ml 5 comes from 



Silbermann &: SmithI (119951 ) 



3. Phase Shift Diagrams and Period Changes 

Period changes of variable stars are usually determined with the use of an 0-C diagram, where 
the O stands for the observed time of a particular phase of the light curve, while the C stands 
for the calculated or predicted time of that same phase. Generally the 0-C curve is calculated 
at a star's maximum phase, and this is the case here as well. Nonlinear 0-C curves indicate a 
changing period, and a linear 0-C relation with a non-zero slope signifies an inaccuracy in the 
adopted period (see, for example. Storm, Carney & Beck 1991). In particular, a second-order 
(parabolic) fit indicates a constant period change rate. From the pulsation equation a change 
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in the period of an RR Lyrae means the structure of the star (specifically, the mean density) 
is changing. Hence, the period changes of RR Lyraes can be connected with different phases of 
horizontal-branch evolution. Although steady evolutionary effects do not explain the large period 
decreases or the random and/or abrupt period changes that are seen in some local and GC RR 
Lyraes, a period increases as a star evolves redward and decreases as it moves blueward. Also, any 
HB star's luminosity evolves monotonically brightward. As an RR Lyrae's luminosity increases so 
does its radius increase and its density decreases. Therefore, a star that evolves both blueward and 
brightward should display a relatively small period change; a star that evolves both redward and 
brightward should display a slightly larger rate of period increase. See §4.2 for more details. 

The period changes of the IC 4499 stars are investigated using the most accurate photometry 
(errors smaller than 0.05 mag). The V and B filters were chosen from among the six available 
{U ,B ,V ,R,I ,D D051) because the V data are the most abundant, and t he early photographic 



plates (needed to extend the baseline) were mostly taken through a B filter. iLee Sz CarnevI (|l999al ) 
show that for the globular cluster M2 the phase shifts between the B and V light curves are quite 
similar: the mean phase shift between the V and B light curves is only 0.0021 (ifc 0.0039 rms error) 



for the 19 RRO stars, and 0.0001 (± 0.0021 rms error) for the 12 RRl stars. lJurcsik et al.1 (|2001 
point out that this same difference is less than 1-2% of the period. 

To construct 0-C diagrams, the V photometric measurements were generally divided into five 
groups of time. The B pho tometric data were d ivided into either four or five groups of time, 
depending on the number of Clement et al. (jl979l ) observations. In general, each group of data 



blocked by time has ~45 photometric values. The data in each time group were folded by the 
period listed in Table [1] and shifted in phase until maxima aligned at = 0. To help fin d where 



0.0, each light curve was fitted with a series of six light curve templates as described in iLavden 



(| 19981 ). The resulting fit was visually examined to ensure that the fitted template accurately 
describes the light curve shapes, especially at maximum light. An example of the quality and 
quantity of data in each time group, as well as an example of the light curve templates' ability to 
align 0=0, are shown in Figure [TJ 

Constructing an 0-C diagram in this manner requires a relatively complete light curve for each 
time group. Hence it can happen that although there may be some data at a given time interval, if 
the data are not evenly spaced about the light curve, no accurate maximum light estimate can be 
made. Although this limits the number of values plotted in the 0-C diagram, it also ensures that 
each point on the 0-C curve diagram is determined very accurately, regardless of whether the RR 
Lyrae light curve shape is slowly changing. For the RR Lyraes in IC 4499 such a procedure works 
well, as the photometry can be easily divided into groups of time with adequate coverage for each 
group (see Figured]). The resulting period change rates were compared to those obtained from one 
of the author's (JMN) SAS program (Nemec, Hazen-Liller &; Hesser 1985) and yielded consistent 
results. 



There are two relatively large gaps in the data, which have the potential to introduce cycle- 
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count ambiguities. Stars with shorter periods and larger period change rates are most prone to this. 
To check for this, each V and B hght curve was divided at HJD = 2,449,000 and a direct period 
was calculated for these two groups using the RR Lyrae period finding approach developed by 
Lavden et al.l (|l999l ). For the stars with the largest period changes and hence those with potential 
cycle-count ambiguities, the difference between the directly calculated periods of the earlier and 
later HJD time groups would indicate whether the period change rates are in the right direction. To 
calculate the 0-C values for these stars, cycle counts were assumed that produced values continuing 
the previous trend in the 0-C diagram. 

The shift required to align the light curve at (j)=0 (phase at maximum) is plotted against the 
mean HJD in Figure [2]l6l A weighted se cond-order (paraboli c) fit as well as the scatter about the 
fit are calculated for each star. Following I Jurcsik et al.l (120011 ) . the stars with a period change rate, 
< AP/At >, of 



I < AP/At > I > 2aAP/At (1) 

are likely to have a monotonic period change. Here P is the assumed period listed in Table [TJ 
Table [U consists of (1) the variable number, (2) the adopted period, (3) the values of /? in d Myr~^, 
(4) the formal error in in d Myr~^, (5) the normalized period variations, a = (3 /Pa in 10~^'^d~^, 
and (6) the ctap/a* x 2 in d Myr~^ for all the variables studied. 

As the above procedure requires a thorough examination of the RR Lyrae light curves, it 
became evident that many stars exhibited Blazhko behavior (amplitude and/or phase modulation, 
see e.g., Jurcsik et al. 2009, Kolenberg 2010). Stars that exhibit the Blazhko behavior must be 
treated with special caution, as disentangling the Blazhko behavior from the period change behavior 
requires a lot of data and is not always straightforward. Although we attempted to solve for the 
period change rates for the Blazhko stars, the errors in (3 are quite large. More data would help beat 
down these errors. We would be extremely grateful to any investigators willing to donate copies 
of proprietary CCD images or calibrated photometric measurements of IC 4499 variable stars to 
enhance our data set for this cluster. The 21 stars (31% of the IC 4499 RRO population) that 
appear to exhibit Blazhko behavior are: VI, V15, V19, V20, V23, V30, V38, V43, V48, V52, V53, 
V57, V58, V61, V7 4, V80, V82, V83, V84, V85, V88. This is a larger percentage than seen in a; 



Cen (20% foun d by I Jurcsik et al 



(47% found bv I Jurcsik et al.ll2009l ^ 



200 ll ). and a smaller percentage than seen in the field RRO stars 



4. Horizontal-Branch Evolution 

4.1. Period Change Rates 

Figure [7] is a histogram of the 31 RRO and the 11 RRl variables that are fitted with parabolas 
that yield period change rates that predict acceptable periods for phasing our data. It is striking 
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that both RRO Lyrae and RRl-type variables tend to have positive period change rates. The 
distribution of the period change rates for the IC 4499 RR Lyrae variables does not scatter evenly 
about zero, but there is a preference for f3 values between 0.0 and 0.1 day Myr~^. There are 22 
RRO stars with positive period change rates and 9 RRO with negative period change rates. For the 
RRl-type variables, only one star with a negative period change rate is seen. The stars not included 
in this figure are the possible RR2("RRe", second-overtone pulsation), RROl ("RRd" double-mode 
pulsation), most of the probable Blazhko variables, stars that have ambiguous O — C curves (V26, 
V32, V38, V96, V106) and stars for which there were not enough data to draw any firm conclusion 
about their long-term period changes (V45, V46). The general shape of the histogram does not 
change upon removing the 12 RR Lyrae variables (6 RRO and 6 RRl variables) with period change 
rates that do not conform to Equation 1. The ratio of stars with increasing periods to total stars 
with calculated period change rates is approximately the same (~70%). 



As seen in other stud ies of period change rates of RR Lyrae stars in GCs {e.g., lReidlll996 



Silbermann Sz Smithlll995l ). there is an excess of large period changes among the RRO stars com- 
pared to the RRl stars. A plausible explanation for larger period changes among fundamental 
pulsators when compared with first overtones is explained in the topology of the instability strip. 
The region where fundamentals are pulsationally unstable is a factor of 3-4 larger (in effective 
temperature) than for the first overtone. This would imply that you have a larger probability to 
detect period changes among fundamental RR Lyrae variables. 



4.2. Period Change Rates From Stellar Evolution 



The period changes discussed above can in principle be used to test models of HB evolution. 
However, the fact remains that many RR Lyrae period change rates have been measured to be an 
order of magnitude larger than predicted from stellar evolution theory in a variety of GCs from 
a variety of sources, and this is the case for the RR Lyrae variables in IC4499 as well. It has 
been shown that when finding period change rates for RR Lyrae stars, there can be "noise" which 
can be an order of magnitude larger and is attributed not to a star's secular evolution but_to, 
e.g., mixing events in the stellar core (jSweigart &: Renzinilll979l ). hydromagnetic events (jStothers 



1980), convection (IStothersll20ldl. or passage throu gh the RR Lyrae instab i lity str ip by stars which 
are pre-ZAHB pulsators (jSilva Aguirre et al.ll2008l ). However. iBono et al.l (Il997bl ) used the metal- 
rich RR Lyrae to show that erratic period changes have timescales significantly shorter than the 
evolutionary effects, so the probability to detect objects during an erratic change is much smaller 
then for the secular change. 

So why is there an abundance of RR Lyrae period change rates seen that are an order of 
magnitude larger than those predicted by HB models? One possible explanation is tha t the p redic- 
tions of the theoretical models are incorrect. This was put forth bv iLe Borgne et al.l (|2007l ) . who 
determined period change rates for local RR Lyrae variables from the GEOS RR Lyrae database-a 
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database containing about 50 000 maximum times from more than 3000 RR Lyrae variables @ Their 
O — C plots were produced from data spanning more than 50 years (and for several objects the 
data exceeded 100 years). A large number of their RRO Lyrae variables showed evidence of linear 
period variations, and the determined period change rates were of the same order of magnitude, 
the median values being j3 = +0.14 d Myr~^ for the 27 stars with increasing periods and (3 = —0.20 
d Myr""*^ for the 21 stars with decreasing periods. They further note that although a /3 value of 
0.30 d Myr~^ is exceptional and the lifetimes of these phases would be very short according to HB 
models (see also Fig ED, in their sample, they had 18 stars with < /3 < 0.30 and 9 stars with (3 > 
0.30. This 2:1 ratio seems too small to be caused by the presence of exceptional cases. 

Another explanation for the abundance of period change rates observed that are an order 
of magnitude larger than predicted by theoretical models, was put forth by iLed ()l99ll ). They 
demonstrate that horizontal branch simulations can reproduce the observed distributions of /3 in 
GCs if random observational errors are superposed on the evolutionary period changes. Hence the 
observed RR Lyrae period change rates can be attributed to evolutionary effects if the random 
observational error is of order ±0.07 d Myr~^ in /?, a value consistent with those suggested by 
the obse rver s . Alt hough the more recent HB models show basically the same behavior as those 
used bv iLed (|l99ll ). recent observationa l investigations present period chang e rates with smaller 



observational errors (Jurcsik et al 



2001 



Nemec 



2004 



Le Borgne et al.l 120071 ). In general, this is 



the case for the observed period change rates presented here. It is difficult to understand where 
and how such ±0.07 d Myr~^ errors in /3 would arise from these newer observations. 

In Figure El the V,B — V CMD for all the stars for which period change rates are determined is 



shown. Th e evolutionary tracks for O.63M0 to O.66M0 from the BaSTI archive (jPietrinferni et al 



2004 



20061 ) are also plotted, where the tracks were calculated with [Fe/H] = —1.62 dex and [a/Fe] 



0.4 dex {e.g., Walker &: Nemec 1996), and a mass loss rate 77=0. 40 The central He exhaust 



ion 



for the adopted chemical composition is represented by the solid por tion of the individual HB 
models. We adopt E{B - V)=0.23, Mv=0.52 and < ^zahb >=17.72 dWalker fc Nemeclll996l ') to 
align the models from the BaSTI archive to the observed CMD. The BaSTI tracks indicate that 
most of the IC4499 RR Lyrae stars have masses of 0.64Mq and 0.65Mq. The general shape of 
the evolutionary tracks changes little assuming a slightly more metal-poor HB model of [Fe/H] 
= —1.82 dex, although the masses that fit the observed colors increases by ~0.03Mq. For both 
of these metallicities, the predicted masses are i n agreement with the p ulsational masses of 0.65- 
O.7OM0 found from the RROl (RRd) variables (jWalker k Nemeclll998l '). However, the [Fe/H] = 
—1.62 dex HB models yield RRO Lyrae periods which are in much better agreement to those found 
observationally for the IC4499 RR Lyrae variables. For the [Fe/H] = —1.82 dex tracks and the 
mass range that fit the observations, the RR Lyrae periods are longer by ~0.1 days. 



^The GEOS database is accessible at |http : //dbrr . ast . obs-mip . f r7] hosted by the Laboratoire d Astrophysique 
de Toulouse- Tarbes, Observatoire Midi- Pyrenees, Toulouse, France. 

The BaSTI archive is available on-line at |http : / /albione . oa-teramo . inaf . it/] 
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in Figure [9] the evolutionary rates of period change in d Myr ^ are plotted as a fun ction of the 



time e lapsed since the ZAHB. To estimate (3, the fundamental pulsation equation of iBono et al 



(11993): 



log Pf = 11.627 + 0.823 log L - 0.582 log M - 3.506 log T^ff (2) 

has been used. The heavy lines represent the observed range of the IC4499 instability strip (0.115 
> {B — V)q > 0.365; see Walker & Nemec 1996) for each track. The evolutionary tracks run first 
redward (see Fig[8|) giving positive period change rates of 0.03 d Myr~^. Soon, blueward evolution 
through the IS occurs, relatively slowly, giving small negative rates of period change. Faster positive 
period changes occur as the redward evolution proceeds, yielding P values of up to ~0.1 d Myr~^. 
These values are obviously much smaller than the observed values for IC 4499. 

In the bottom panel of Figure [10] the histogram of the evolutionary period change rates are 
presented from the [Fe/H] = —1.62 dex HB tracks discussed above, assuming 15% of the RR Lyrae 
variables have O.63M0 and 0.65Mo and 35% of the variables have O.64M0 and O.65M0 (see FigE]). 
To simulate observational errors, the evolutionary period change rates are enhanced by an order 
of magnitude (~ib0.07 d Myr~^). Because the baseline of the IC4499 observations is quite long, it 
is rather straightforward to observe if the period of an RR Lyrae star either increases or decreases 



over time. For this cluster, the addition of random observational errors, as introduced by iLee 



(jl99ll ). would not make sense. Instead we introduce a systematic error, as could be imagined for 



example, if the errors on the different 0-C points were underestimated and h ence the weighted 



parabolic fit would yield greater coefficients than realistic. As with iLed (|l99lh . we find that the 
observed period change rate histogram is very similar to that predicted by HB models. Repeating 
this procedure using the [Fe/H] = —1.82 dex HB tracks with masses appropriate for the observed 
IC4499 observations, a similar shaped histogram is obtained but with a slightly wider base. This is 
not surprising given the similarity in the shape of the evolutionary tracks at these two metallicities 
and that the evolutionary period change rates of more metal poor RR Lyrae variables tend to be 
larger. 

In the top panel of Figure [TOl the period change rates of the IC 4499 RR Lyrae are compared to 
those in M3, because M3 has a similar HB parameter and metallicity to IC 4499. The histograms 
look very simil ar. In particular, the p eriod change rates are of the same order of magnitude. As 



pointed out by iLe Borgne et al.l (120071 ) since almost all of RR Lyrae period change rates have the 
same order of magnitude they should have a common origin. It would be very unlikely that all 
the 42 RR Lyrae variable in IC 4499 and 37 RR Lyrae variables in M3 for which period variations 
are determined, are undergoing period variations caused by various instabilities, tidal effects, etc. 
However, we cannot rule out the possibility that some individual stars in the 42 IC 4499 and 37 
M3 RR Lyrae with period change rates are merely in a particular transition phase. It is much 
more likely that the period changes for the RR Lyrae i n IC 4499 and M 3 are caused by a more 



general phenomenon, with the immediate answer from iLe Borgne et al.l (|2007l ) being long-term 



stellar evolution. Our period change rate values are similar to those reported by (the more recent 



- 11 - 



papers only) lJurcsik et al.l (120011') in the case of u Cen R R Lyrae stars, iLe Borgne et al. (2007) in 
the case of the field RR Lyrae, ICorwin &: CarnevI (|200ll ) in the case of M3, and iNemed (|2004l ) in 
case of NGC 5053, and discussed as evolutionary changes. 

In conclusion, we can not prove that the period changes are due to evolution since the com- 
parison with the models are poor. The HB models predict period change rates that are an order of 
magnitude sma l ler th an what is observed. Given the distribution of the errors in the analysis, the 
analysis bv lLed ( 199 ll ). which suggests an observational error of ±0.07 d Myr"'^ in /3, does not seem 
valid, either. Recent observed period changes have errors smaller than what is required to match 
the observations and models in a satisfactory manner. The theoretical models should be fine tuned 
to match the observed rates in a better way, now that there are extensive databases available that 
allow period change determinations. 



4.3. The Period- Amplitude plane 

The Period- y- Amplitude (PAy) plane is often used to distinguish between Ool and OoII 
globular clusters. Here we explore how the period change rates of the RR Lyrae variables in Ool 
and OoII globular clusters compare in the PAy plane on a star-by-star basis. 

Figure [TT] shows the PAy plane for the IC4499 RRO variables. As IC4499 is an Ool cluster, 
more RR Lyrae variables fall on the Ool line. If RR Lyraes in Ool globular clusters were evolving 
from red to blue, their periods would be predominantly decreasing, and vice versa for OoII GCs. 
However, most of t hese Ool RR Lyraes have positive period change rates. Here the Ool and OoII 
lines are defined by Clement &: Shelton (jl999l ): the Ool line is the least-squares fit to RRO stars in 



the Ool prototype GC M3, and the OoII line is the least-squares fit to RRO stars in the OoII GC 
M9. 

Blazhko variables tend to introduce scatter in the PAy plane due to their non-repeating light 
curves, and hence their amplitudes vary over timescales longer than the basic pulsation period. The 
Blazhko variables are denoted by triangles, and the RRO variables for which period change rates 
could not be determined are also indicated. By and large, these stars fall in the same region of the 
PAy plane as the stars with well determined period change rates. Hence a bias in which stars were 
sampled is minimal. 

The PAy diagram of the RRO variables in the O osterhoff I GCs M5, M3 and NGC 70 06 is 



shown in Figure [T2l The period change rates come from lReidI ()1996l ). ICorwin &: Carney! (120011 ) and 
Wehlau, Slawson & Nemec (1999). Again, there is a preponderant number of stars evolving from 
blue to red; on average, there are 1.5 times as many RRO variables with positive period change 
rates compared to negative period change rates. There does not appear to be a trend of where the 
RRO variables with either positive or negative period change rates lie. 



Stars evolving from red to blue (negative /3 values) along the horizontal branch may have a 
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period change rate app roximately four times slower than those moving in the reverse direction 
( Lee &: Demarquelll99d ). Hence, it could be that many stars with no measured period change rates 



have negative (3 values that are too small to be measured. Although this possibility can not be 
ruled out, even if all the stars with no determined period change rates had a negative period change 
rate, it is still striking how many RR Lyrae variables with positive period change rates exist in 
the Ool clusters. HB stellar models with metallicity suitable for Ool GCs show the blue excursion 
responsible for the negative period change only for a short interval of their HB lifetime. As seen in 
§4.2 for the IC 4499 variables, the blue excursion is ~5 Myr (6%) of their HB lifetime. 



In Figure HSl the RRO variables in the OoH GCs M2, NGC 5466, M15, M22 and NGC 5053 
are shown on a PAy diagr am. The period change rate s come fromlLee &: Carneyl (Il999al ). Corwin, 
Carne y & Nifong (1999), Isilbermann Smithl lj 19951 ). Iwehlau fc Sawer Hog J ( 19781 ) and iNemec 
(|2004l ). respectively. Although there are fewer RRO variables with period change rate determina- 
tions, it is clear that there are proportionately more stars with positive period change rates than 
negative ones. As with the Ool globular clusters, most RRO variables appear to have positive period 
change rates. If the period change rates can be attributed to evolution, this would mean both Ool 
and OoII RR Lyrae variables would be evolving from blue to red. 

Recently, high quality period change rates based on a plethora of data and long (~100 year) 
baselines hav e been dete rmined for RR Lyraes in the globular cluster uj Cen and in the field by 
Jurcsik et al.l (120011 ) and lLe Borgne et al.l (|2007l ) . respectively. Neither the Milky Way field, nor cj 
Cen is formally seen as an Ool- or OoII-type object. Rather, it has been suggested that oj Cen may 
be th e remnant of a dwarf galaxy that was accreted by the main bod y of the Milky Way long ago 



{e.g., iNorris et al. 



1996 



Metz et al 



2005 



Da Costa fc Colemanll2008l ). The field R R Lyraes in the 



Gala. c tic halo seem to di s play the separation into the two different Oosterhoff g roups (jSuntzeff et al 



1991 



Bono et al.lll997cl : iMiceh et al.ll2008l : Ide Lee 



2008 



Szczvdel et al.ll2009l l. 



The PAy stars from these two studies are shown in Figure [Ml Among the field RR Lyraes, 
blueward evolution (decreasing periods) is as common as redward (increasing periods). For uj Cen, 
period increases dominate the RR Lyrae population. The RRO variables with negative period 
change rates again do not fall preferentially on either the Ool or OoII line, but are scattered 
randomly in the PAy plane. Again, there does not seem to be a relation between position on the 
PAy plane (Oosterhoff type) and period change rate (evolution?). 

To quantify any differences between period change rates and Oosterhoff types, the mean, 
median and mode period changes for stars in two difi^erent parts of the PAy diagrams are derived. 
We define the Ool section, Ool — PAy, as: 

Av < (-0.73/0.13) P + 4.3 (3) 

which consists of variables closer to the Ool line in the PAy diagram. The OoII section, OoII— PAy: 

Av > (-0.73/0.13) P + 4.3 (4) 



consists of the variables closer to the OoII line. 
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Table 2 displays the results for the clusters studied above and the Oosterhoff sections with 
~10 RRO Lyrae variables. In Table 2, we list (1) the particular population, GC or group of GCs, 
(2) the Oosterhoff section as defined by Equations 3 and 4, (3) the mean period change rate and 
the error in the mean, in d Myr"-*^, (4) the median value of /3 in d Myr"-"^, (5) the mode of /3 in 
d Myr~^, where the bin size used in determining the mode is typically 0.05 d Myr~^, and (6) the 
number stars in each Oosterhoff section. In general, the RRO Lyrae variables in all the populations, 
GCs and Oosterhoff sections tend to have increasing period change rates. Moreover, the size of the 
period change rates are very similar for all the different Oosterhoff sections. There is no indication 
of correlation with period change rate and position in the PAy plane. 



De-reddened Colors 



A search for differences in the intrinsic temperatures between Ool and OoII RR Lyraes is also 
carried out using their ( V-R) colors at minimum lig ht, (V-R)m.in,o- Not only is (V-R) known 



to be a good tempera ture indicator for variabl e stars (jBarnes &: Evans 



1976 



Barnes et al 



1977 



Manduca Bel]|ll98ll ). but lKunder et al.l (|20ld ) also showed that the ( V~R)min,o for RRO variables 
is constant to 0.02 mag and is largely independent of period, amplitude and metallicity. Therefore, 
( V-R)min,o can be used as in indicator to probe for differences in the temperatures of Ool and 
OoII variables. For example, it is well-kno wn that Ool clusters tend to be more metal-rich and 
host fainter RR Lyraes than OoII clusters (|Caputdl200d ). As the metallicity has an effect on the 
absolute magnitude of an RR Lyrae, it has been difficult to disentangle whether the metallicity 
difference alone is affecting the brightness differences, or whether there are differences in the intrinsic 
magnitudes of RR Lyrae variables in Ool and OoII globular clusters caused by something other 
than just metallicity (like evolution or helium). Unfortunately, there is little V- and /2-band data 
for RR Lyrae variables in globular clusters to find their minimum light ( V-R) color. Using new V 
and R data of RR Lyrae variables in IC 4499, we investigate ( V-R)min,o for the Ool cluster IC 4499 
compared to the OoII cluster M15. 



Based on main-sequence fitting and a period-shift analvsis. iLee fc: CarnevI (jl999bl ) showed that 
the RR Lyrae variables in the OoII cluster M2 are ~ 0.2 mag brighter than those in the Ool cluster 
M3, despite the fact that these clusters have similar [Fe/H]. If their results can be generalized to all 
OoII and Ool RR Lyrae stars, this would mean that OoII stars have brighter absolute magnitudes 
than Ool stars, even at the same [Fe/H] . 

If the luminosities of RR Lyrae stars depend on Oosterhoff type, are there other intrinsic 
differences between the RR Lyraes in Ool and OoII type clusters? More importantly, could such 
intrinsic differences shed light on the origin of the Oosterhoff groups? This is addressed by inves- 
tigating the dereddened minimum light (V-R), {V-R)minfli of the fundamental mode RR Lyrae 
stars in the Ool cluster IC 4499 and the OoII cluster M15. Because color is a te mperature indicator , 
and because {V-R)min,o is independent of period, amplitude and metallicity (IKunder et al.ll20ld ). 
potential temperature differences inherent in the different Oosterhoff groups would be elucidated. 
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Figure [T5] shows the dereddened minimum light color of the IC4499 and M15 RRO variables 
as a function of radial distance from the center of the cluster. The radial distance is used to rec- 
ognize possib le differential reddening effects. The adopted color excess for M15 is E(B — V)=0.10 



Harris 



19961) ■ which is a lso supported by analysis of deep Hubble Space Telescop e photometr y 
Recio-Blanco et al.l[2005l ). The adopted color ex cess for IC 4499 is E(B - V)=0.23 (|Harrislll996l '). 
This value is between E(B — V)=0.22, found by IWalker &: Nemed (119961) using four different es- 
timates of reddening to the cluster, and E(B — V)=0.24, found by (IStormI l2004l ) using the mean 



V- and IT-band photometry of RR Lyrae variables in IC 4499. The minimum-light color is found 
by averaging the ( V-R) color over the phase interval 0.5-0.8, and the ( V-R)min colors for the 
IC4 499 are presented in Table 3 . For the M15 RR Lyrae variables, only those with light curves 
that ISilbermann fc SmithI (jl995l ) consider "good" were used, and these stars have photometric 
standard errors of ~0.015 mag in ( V-R). To convert from E(B — V) to E(V — R), it is assumed 
that E(V - R) = 0.77 E(B - V) from Cardelh, Clayton & Mathis (1989). 

The average dereddened (V-R) at minimum light, {V-R)min,o for the M15 and IC4499 RRO 
variables is 0.28 it 0.02 and 0.27 it 0.02, respectively, where 0.02 is the dispersion about the 
mean. The standard error of the mean is 0.003 and 0.007 respectively, although these numbers 
do not include any possible systematic errors in the photometric calibration of either study, which 
could — in our experience — be of order 0.01 mag. Not only are these values almost identical to the 
{V-R)min,o for local RR Lyrae stars (jKunder et al.ll20ld ). but these values are within the errors 
of each other. Any differential reddening is ncgligiblG compcircd. to the errors in <C ^~R^min,0' 

At 

least for this Ool and OoII cluster pair, no difference in the dereddened colors at minimum light is 
seen in RR Lyrae variables at the 0.02 mag level. 

This result suggests that the intrinsic color differences between the Ool and OoII RR Lyraes 
are either very small, or do not occur at minimum light. As ( V-R) color is a good indicator 
of temperature {e.g., Barnes & Evans 1976), our result also suggests that at minimum light Ool 
and OoII RR Lyraes have very similar temperatures. A 0.02 mag uncertainty translates into a 
temperature differei ice between the Ool and OoII RR Lyraes of ~ 100 K at minimum light. This is 
calculated from the lBarnes et al.l (|l977l ) relationship relating effective temperature to dereddened 
Johnson ( V-R) from Cepheid pulsation, and agrees with model atmosphere parameters appropriate 
to RR Lyrae stars from Manduca & Bell 1981 (their Table 1) for both Johnson and Cousins color^. 

In order to fully understand observed prope rties of RR Lyrae variables , it is necessary to study 
those properties as a function of phase (see e.q.. lKanbur &: Fernanddl2005l ). Although the shape of 
the light curve affects the estimate of the mean magnitude and colors (Bono, Caputo &: Stellingwef 
1994) it appears that at minimum light, RR Lyraes are remarkably similar. It would be worthwhile 
to perform a more comprehensive analysis of the various period-color and amplitude-color relations 
as functions of pulsation phase and Oosterhoff type. This is beyond the scope of this paper, as 
detailed knowledge of how the metallicity, period and amplitude affect other pulsation phases would 



*We note that the Johnson R used in 



Barnes et al. 



( 19771 ) is not the same as the Kron Cousins R used here. 
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be needed. 

Given the homogenous photometry and reddening of the RR Lyraes in IC4499, it is promising 
to search for any ( V-R)m.in,o dependence on period or amphtude. A weighted least squares analysis 
of ( V~R)m.in,o as a function of period reveals a slope significant at the 2 a level (0.09 it 0.04 mag 
day"^). Although not formally significant, this may suggest ( V-R)min,o does have some dependence 
on period at the ~0.01 mag level. A weighted least squares analysis of ( V-R)min,o as a function of 
amplitude reveals no correlation (the slope found is —0.008 ± 0.045 mag mag~^). 



6. Discussion 



A very basic prediction of stellar evolution and pulsation theory is that the period of a star 
changes with time. The luminosity and temperatu re of a s t ar wi ll slowly change, and correspond- 
ingly its mean density will also change. From the iRitter I (|l879l ) period-mean density relation of 
pulsation theory, P \/< P > ~ const., the period should change inversely with the density. The 
models of lLed (Il99ll ) clearly show a predicted trend of increasing period change rates for bluer HB 
types (which are generally OoII clusters) and decreasing period chang e rates for redd er H B types 
(gener ally Ool clusters). Among the more recent papers, we mention ICatelanI (j2009l ) and iNemec 
(j2004l ) who fit these models to the observed GC RR Lyrae period change rates. Although in general 
the large uncertainties in the mean period change rates do overlap these models, the chi squared 
values are not consistent with a good fit (x^~0.08). The above analysis indicates why this may be 
the case-there is a predisposition for positive period change rates in both Ool and OoII RR Lyrae 
variables in the Milky Way GCs. No cluster of either Oosterhoff group exhibits an abundance of 
RR Lyrae variables with decreasing periods. 

The RRO Lyrae variables in u Cen and in the field confirm this result. Although most RR 
Lyrae stars in u Cen are known to be similar to those in OoII clusters, a smaller population are 
similar to RR Lyraes in Ool clusters (Butler, Dickens, & Epps 1978; Caputo 1981; Clement &: 
Rowe 2000). These two populations are easily distinguished in the PA plane; however, there is 
no evidence of a correlation of period change rate as a function of their position on the period- 
amplitude diagram. The same is true for the field RR Lyrae variable s , which are also shown to be 



a mix of RR Lyraes wit h Ool and OoII properties (IMiceli et al.ll2008l : Ide Ledl2008l : ISzczygiel et al 



20091 : iBono et al.lll997d ). Assuming that the period-amplitude diagram can be effectively used to 
classify RR Lyrae stars into an Oosterhoff type, this means that evolution and Oosterhoff type 
are not correlated. One of the most promising ii iterpretations of the double ma.in seq uence in the 
CMD of uj Cen is a helium abundance variation (|Norris et al.ll2004l : iPiotto et al.ll2005l ). which may 
suggest that helium abundance va riations contributes to the Oosterhoff types. Prom high-resolution 
spectroscopy, ISollima et all (|2006l ) found that the helium enrichment seen in the oj Cen blue main 
sequence stars is absent in the metal intermediate RR Lyrae population. Hence at least for uj Cen, 
two populations with similar metallicities but very different helium abundances seem to coexist 
within the cluster. Further, as shown from nonlinear pulsation models including a non-local time- 



- 16 - 



dependent treatment of convection from e.g.. iMarconil (|2009l ) and lBono et al.l (|l997l ). differences in 
helium can significantly contribute to the position of an RR Lyrae in the period-amplitude diagram. 

If differences in HB evolution are the explanation for the different Oosterhoff types, RR Lyrae 
stars in OoII cluster would begin their lives at higher temperatures and evolve to lower temper- 
atures, exhibiting positive period change rates as they evolve. From the minimum light colors 
of the RR Lyrae variables in the Ool cluster IC4499 and the OoII cluster M15, no difference in 
temperature at minimum light is seen to within ~ 100 K. It would be interesting to investigate any 
temperature and color differences between RR Lyrae variables in different Oosterhoff type GCs as 
a function of other pulsational phases. 

Recently two metal-rich GCs, NGC 6388 and NGC 6441 ([Fe/H] 0.6 dex), were found to host 

a number of RR Lyrae variables with properties th at do not follow the period-metallicity trends of 
the other two Oosterhoff groups (jPritzl et allboOlh . The RR Lyrae variables in these GCs occupy 
the OoII region in the PA plane and have mean periods similar to OoII GCs, but shifted to even 
longer periods. However, they have metallicit ies similar to, or ev en higher than, Ool GCs. Termed 
as "Oosterhoff III" GCs, it has been shown bv lPritzl et al.l (|200d ) that the RR Lyraes found in these 
GCs can not be totally accounted for by evolutionary effects. Although these GCs are peculiar in 
their own right, this result is consistent wi t h our findings that evolution and Oosterhoff type are 
not correlated. As detailed by Busso et al. (j2007l ). a possible solution to the peculiar morphology 
of the HB in both clusters is the p resence of multiple stellar populations with different initial He 
contents. ICaloi D'Antonal (|2007l ) have also shown that the mean period of RR Lyrae in these 
two strange clusters can be reproduced by assuming a strong He enhancement. 



It is also worth mentioning NGC 4147, an Ool globular cluster (IStetson et al.l 120051'). This GC 
has one of the bluest HBs and is one of the most metal-poor Ool clusters. IContreras et al.l (j2005l ) 
point out that an Ool classification for NGC 4147 is inconsistent with the theoretical paradigm 
that RR Lyrae stars in Ool globular clusters are relatively unevolved objects, whereas those in 
OoII globulars are evolved from a position on the blue ZAHB. Further, this globular cluster, as 
well as the Ool GCs NGC 6171 and M62 have a large number fraction of RRl variables to RRO 
variables, which is in stark contrast to the theory that Ool GCs are less evolved than OoII GCs 
(jContreras et al.ll2005l ). 



Seventeen doub le-mo de RR Lyrae sta r s (RR Ol stars) in IC4499 have been identified by 
Clement et al. (jl986l ) and Walker &: Nemec ( 1996 ). Because the period ratios of RROl stars have 
a strong dependence on mass, derivation of masses for R R Lyrae stars c an be obtained {e.g., 
Cox, King & Hodson 1980, Cox, Hodson & Clancy 1983). Iciement etall (Il986l l found that the 
IC4499 RROl stars are similar to those stars in the Ool cluster M3, but considerably different 
from RROl stars found in various OoII systems. In the period ratio diagram (used for stars that 
pulsate with at least two periods simultaneously; the ratio of the periods is plotted as a func- 
tion of long er period), the RROl s tars divide into two groups, split according to Oosterhoff type. 
Since then, (jClementini et al.ll2004l ) have found two M3 RROl stars with period ratios so low, that 
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they lie well separated from all previously known RROl stars in the Petersen diagram. This large 
spread in period r atio o f the M3 double-mode pulsators is most likely due to spread in mass, and 



Clementini et al.l (|2004l ) suggest this anomalous spread in mass could be caused by mass-transfer 



in binary systems, helium enhance ment, or, less likely, varying a-element enhancement among the 
M3 stars. IClementini et al.l (|2004l ) also point out that the scenario that Ool clusters such as M3 
are evolving blueward does not fit well, since both blueward and redward evolution seems to occur 
among the M3 double-mode RR Lyrae stars. 



Conclusions 



Period change rates for the RR Lyraes in IC 4499 are calculated and it is found that there is a 
preference for the RR Lyrae variables in IC 4499 to have positive period change rates. These period 
change rates are an order of magnitude larger than predicted by HB evolution models. Although a 
systematic error of ±0.07 d Myr~^ in /3 would reconcile difference between theory and the period 
change observations presented here, the error analysis in the period change rate determination does 
not easily lend itself t o such a larger error. Recent period chan ge rate determinations from time 
baselines of 100 years (jLe Borgne et al.l 120071 : 1 Jurcsik et al.ll200ll ) equally present challenges to the 
period change rates predicted by the HB models. 

If the period change rates can be attributed to stellar evolution, this surplus of period increases 
is in contrast to the explanation that the RR Lyraes in Ool clusters are less evolved than those in 
Ool clusters and preferentially traverse the instability strip from red to blue. Comparisons with the 
Ool clusters NGC 7006, M3 and M5 also show a preponderance for RR Lyraes with positive period 
change rates, which may imply that the RR Lyraes in the Ool clusters also do not preferentially 
evolve from red to blue along the horizontal branch. Comparisons of the period change rates of 
the OoII clusters M15, M2 and NGC 5053, the field RRO variables, and the mixed cluster oj Cen, 
imply this as well. Although OoII RRO variables cluster along a specific plane in the PAy diagram, 
there is no clear preference for these RR Lyrae variables to have either positive or negative period 
change rates. Hence if the period change rates are an indication of evolution, differences in the 
evolutionary states of the Ool and OoII RR Lyrae variables are not seen. This would strongly 
suggest that age is not the primary answer in explaining the Oosterhoff types. 

We have discussed the dereddened ( V-R) colors and temperatures at minimum light of the 
Ool cluster IC4499 and the OoII cluster M15. It is found that to within 0.02 mag, there is 
no difference in the minimum light dereddened colors of the RRO variables from these two clus- 
ters. This is an especially useful property given that there are studies showing that at similar 
[Fe/H] metall icities, the absolute niagnitudes of RR Lyrae star s in Ool and OoII clusters can differ 
by ~0.2 mag (|Lee &: Carnevlll999bl : iKunder &: Chaboveiil2009l ). In terms of Tgff, this suggests that 
the temperatures of the RR Lyraes in Oosterhoff I and Oosterhoff II globular clusters are within 
~100 K at minimum light. If variables in OoII clusters are more evolved than those in Ool clusters, 
this temperature similarity must be explained by a difference in mass, luminosity or Helium content 
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among possible effects. 
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Table 1:: Period Change Rate of IC 4499 RR Lyrae Stars 



Number 


^^^^ J\ / J\ 

Ir'erioa(a) 


p (d JVlyr ) 




ct (1/lU dj 


2 X o-AP/At(d Myr ) 


2 


0.4936367 


0.12 


0.03 


6.49 


0.12 


3 


0.4832464 


0.05 


0.02 


2.61 


0.07 


4 


0.6236188 


0.18 


0.05 


7.77 


0.36 


5 


0.5568808 


-0.30 


0.08 


-14.90 


0.04 


6 


0.5778506 


8.7 


0.1 


238.3 


1.1 


9 


0.7096094 


0.01 


0.02 


0.016 


0.02 


11 


0.6314811 


0.04 


0.03 


1.65 


0.01 


12 


0.5947952 


0.17 


0.06 


8.01 


0.03 


13 


0.5115214 


-0.29 


0.05 


-15.52 


0.36 


16 


0.5597099 


0.10 


0.07 


4.65 


0.03 


17 


0.4980768 


0.38 


0.11 


20.83 


0.22 


24 


0.5165604 


-0.31 


0.08 


-16.17 


0.07 


25 


0.6021555 


-0.03 


0.11 


-3.76 


0.06 


27 


0.5067589 


0.52 


0.02 


28.310 


0.09 


28 


0.5827866 


0.21 


0.10 


9.866 


0.09 


29 


0.3626859 


0.03 


0.04 


2.567 


0.08 


34 


0.4935571 


0.08 


0.133 


4.55 


0.06 


37 


0.5648090 


0.33 


0.03 


16.19 


0.06 


44 


0.5253868 


-0.48 


0.06 


-24.96 


0.03 


46 


0.5747031 


-3.5 


0.4 


-95.9 


0.13 


47 


0.5832672 


2.1 


0.1 


57.3 


0.8 


49 


0.4991130 


1.1 


0.1 


29.9 


0.7 


50 


0.5657635 


1.59 


0.02 


43.5 


0.01 


51 


0.3551330 


0.56 


0.16 


42.94 


1.61 


54 


0.7364995 


-0.91 


0.22 


-33.6 


0.02 


55 


0.3585305 


0.16 


0.03 


12.04 


0.92 


64 


0.6361433 


-0.06 


0.16 


-2.71 


0.07 


66 


0.6260572 


0.22 


0.03 


9.62 


0.10 


69 


0.340198 


4.25 


0.51 


116.3 


2.27 


70 


0.5370717 


1.9 


0.05 


53.1 


0.32 


72 


0.6813507 


0.32 


0.09 


12.90 


0.02 


74 


0.6023051 


0.95 


0.07 


43.18 


0.04 


76 


0.6416066 


0.84 


0.21 


35.93 


0.02 


77 


0.3546419 


-0.20 


0.01 


-15.52 


0.09 


85 


0.5939029 


0.46 


0.32 


-20.56 


0.07 


89 


0.3553893 


0.12 


0.01 


9.24 


0.12 


92 


0.3547447 


0.17 


0.02 


13.35 


0.25 


95 


0.3626458 


0.14 


0.02 


10.27 


0.23 


98 


0.3556421 


0.51 


0.06 


39.47 


0.02 


103 


0.3529834 


0.11 


0.03 


8.76 


1.05 


108 


0.6414954 


-0.66 


0.16 


-27.96 


0.12 


111 


0.3573436 


0.14 


0.01 


11.03 


0.02 
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Table 2: Period Change Rates as a function of position on the PAy diagram 



Cluster 


Oosterhoff section 


mean (d Myr median (d Myr ^) 


mode (d Myr -"-) 




NGC 7006 


Ool- 


-PAv 


0.04 ± 0.02 


0.03 


0.03 


34 


M3 


Ool- 


-PAv 


0.05 ± 0.04 


0.11 


0.11 


33 


M15 


Ool- 


-PAv 


-0.03 ± 0.04 


0.02 


0.01 


13 


IC 4499 


Ool- 


-PAv 


0.02 ± 0.20 


0.17 


0.08 


19 


Ool GCs 


Ool- 


-PAv 


0.04 ± 0.02 


0.03 


0.03 


80 


OoII GCs 


Ool- 


-PAv 


-0.08 ± 0.11 


0.02 


0.05 


7 


ui Cen 


Ool- 


-PAv 


0.05 ± 0.10 


0.03 


0.07 


13 


Field 


Ool- 


-PAv 


-0.04 ± 0.05 


-0.07 


-0.07 


24 


Ool GCs 


OoII- 


-PAv 


-0.09 ± 0.15 


0.04 


0.2 


10 


OoII GCs 


OoII- 


-PAv 


0.11 ± 0.04 


0.08 


0.15 


27 


oj Ccn 


OoII- 


-PAv 


-0.26 ± 0.48 


0.13 


0.18 


42 


Field 


OoII- 


-PAv 


0.02 ± 0.11 


0.08 


0.07 


10 



Table 3: Minimum Light Colors the IC 4499 RRO Lyrae Stars 



Star Number {V — 




1 


0, 


.46 


± 


0, 


.01 


2 


0, 


.46 


± 


0, 


.02 


3 


0, 


.45 


± 


0, 


.02 


4 


0, 


.47 


± 


0, 


.01 


5 


0, 


.46 


± 


0, 


.01 


6 


0, 


.47 


± 


0, 


.01 


7 


0, 


.46 


± 


0, 


.03 


9 


0, 


.45 


± 


0, 


.02 


11 


0, 


.46 


± 


0, 


.01 


12 


0, 


.46 


± 


0. 


.03 


13 


0, 


.46 


± 


0, 


.02 


14 


0, 


.46 


± 


0, 


.01 


15 


0, 


.46 


± 


0. 


.01 


16 


0, 


.44 


± 


0. 


.01 


17 


0, 


.46 


± 


0, 


.01 


19 


0, 


.48 


± 


0. 


.06 


20 


0, 


.46 


± 


0, 


.01 


23 


0, 


.43 


± 


0, 


.02 


24 


0, 


.47 


± 


0. 


.03 


25 


0, 


.48 


± 


0. 


.01 


26 


0, 


.49 


± 


0, 


.01 


27 


0, 


.43 


± 


0, 


.01 


28 


0, 


.45 


± 


0. 


.01 


30 


0, 


.42 


± 


0, 


.05 


34 


0, 


.41 


± 


0, 


.05 


36 


0, 


.39 


± 


0. 


.07 


37 


0, 


.47 


± 


0. 


.03 


38 


0, 


.44 


± 


0, 


.01 


44 


0, 


.42 


± 


0, 


.02 


47 


0, 


.49 


± 


0. 


.01 


48 


0, 


.42 


± 


0, 


.02 


49 


0, 


.41 


± 


0. 


.04 


50 


0, 


.42 


± 


0. 


.03 


52 


0, 


.44 


± 


0, 


.02 


41 


0, 


.45 


± 


0, 


.02 


57 


0, 


.47 


± 


0, 


.02 


58 


0, 


.42 


± 


0, 


.02 


61 


0, 


.46 


± 


0, 


.02 


64 


0, 


.46 


± 


0. 


.02 


70 


0, 


.48 


± 


0, 


.03 


72 


0, 


.46 


± 


0, 


.01 


74 


0, 


.45 


± 


0. 


.01 


82 


0, 


.42 


± 


0, 


.03 


83 


0, 


.46 


± 


0, 


.01 


84 


0, 


.47 


± 


0. 


.02 


88 


0, 


.47 


± 


0. 


.01 


106 


0, 


.47 


± 


0, 


.01 


108 


0, 


.44 


± 


0, 


.02 


112 


0, 


.48 


± 


0. 


.02 
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Fig. 1. — An example of the construction on a typical 0-C diagram using the F-band data of RR 
Lyrae V9 in IC4499. The five different light curves show the IZ-band data (represented by blue 
dots) in the five different time divisions. The RR Lyrae template used to find 0=0 is over-plotted 
in red for each light curve. The HJD range is indicated above the lightcurves. The light curve with 
the most recent photometry, shown in the upper right hand corner, consists of data from only two 
nights. The dashed line shows the weighted 2nd order polynomial fit. 
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Fig. 2. — The 0-C diagrams for the RR Lyrae variables in IC 4499. The triangles (blue) represent 
the phase shifts derived from the S-band data, and the open circles (red) indicate the phase shifts 
derived from the F-band data. 
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Fig. 3. — The 0-C diagrams for the RR Lyrae variables in IC4499. Symbols are the same as in 
Figure El 
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Fig. 4. — The 0-C diagrams for the RR Lyrae variables in IC 4499. 
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Fig. 5. — The 0-C diagrams for the RR Lyrae variables in IC 4499. 
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Fig. 6. — The 0-C diagrams for the RR Lyrae variables in IC4499 with relatively large period 
change rates. Note an increase in the scale of the vertical axis to encompass these period change 
rates. 
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Fig. 7. — The period change rate of 31 RRO and 11 RRl-type variables in IC4499. The bin size is 
0.4 d Myr~^. The lower panel presents the same data as the upper panel, but zooming in around 
the peak of the distribution. Here the bin size is 0.15 d Myr~^. 
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Fig. 8. — V,B — V CMD of the RR Lyrae for which period change rates were determined. The RRl 
variables are indicated as triangles and the RRO variables are indicated by circles. The stars with 
positive period change rates are indicated by plus signs, the stars with negative period change rates 
are indicated by minus signs. Stars with period change rates that do not conform to Equation 1 
are indicated by light plus and minus signs. BaSTI HB evolutionary model models are plotted, 
with [Fe/H]=— 1.62 dex, alpha-enhanced. See text for details. 
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Fig. 9. — The evolutionary period change rates (/3 in d Myr~^) predicted for the horizontal branch 
stars belonging to IC4499 as a function of the time elapsed since the ZAHB. The heavy lines 
indicate the approximate range of the instability strip in IC4499, as observed bv lWalker &: Nemec 
(jl996l ). See text for details. 
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Fig. 10. — Top: Comparison of the observed period change rates of IC4499 presented here and 



M3 from lCorwin &: Carney! (j200ll ). Bottom: Comparison of the observed period change rates with 
those predicted from evolutionary models, where the evolutionary models include an observational 
error of ± 0.07 d Myr-^ The histo grams have been normalized. 
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Fig. 11. — Period-amplitude relation for the RRO variables in IC4499. The stars with positive 
period change rates are indicated by plus signs, the stars with negative period change rates are 
indicated by minus signs, and the RRO Lyrae variables for which no period change rates were 
determined are circles. Stars with period change rates that do not conform to Equation 1 are 
indicated by light plus and minus signs. The triangles denote stars that show evidence for the 
Blazhko effect. The least-squares fit to RRO st ars in the Ool prototype GC M3 and the least- 
squares fit to RRO stars in the OoII GC M9 from IClement &: SheltonI (jl999l ) are indicated by solid 
lines. 
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Fig. 12. — Period-amplitude relation for the RRO variables in the Ool globular clusters M5, M3, 
and NGC 7006. The stars with positive period change rates are indicated by plus signs, and the 
RRO variables with negative period change rates are indicated by minus signs. The stars for which 
the period changes are un known (or very unc e rtain ) are indicated by open circles. The Ool and 
Ooll lines are taken from IClement &: SheltonI (jl999l ). determined from RRO stars in the Ool GC 
M3 and the Ooll GC M9. 
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Fig. 13. — Period-amplitude relation for the RRO variables in the OoII globular clusters M15, M2, 
and NGC 5053. The symbols are the same as in Figure [T2l 
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Fig. 14. — Period-amplitude relation for the RRO variables in the mixed population u Cen and the 
field RRO variables. The symbols are the same as in Figure [T2l 
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Fig. 15. — The dereddened minimum light ( V-R) color of the RR Lyrae variables in the Ool 
cluster IC 4499 and OoII cluster M15 is plotted against r, the distance from the cluster center. The 
mean {V-R)min,o for both clusters is indicated by a dashed hne and is essentially the same as the 
( V~R)min,o of local field RR Lyrae stars. 



